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We present a protocol to simulate the quantum correlation implied by non maximally entangled 
two qubit states, in the worst case scenario. This protocol makes a single use of PR-box and a single 
use of Millionaire box (M-box). To the best of our knowledge, the resources used in this protocol 
are weaker than those used in previous protocols and are minimal in the worst case scenario. 
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One of the most intriguing features of quantum physics is the non-locality of corre- 
lations obtained by measuring entangled quantum particles. These correlations are 
I nonlocal because they are neither caused by an exchange of a signal, as any hypo- 
thetical signal should travel faster than light, nor are they due to any pre-determined 
agreement (shared randomness) as they break Bell's inequalities flf. A natural way 
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^ to understand these correlations is to classically simulate them ( known as simula- 
1 — '! tion of entanglement) using minimal resources. Obviously, this cannot be done using 
only local resources, that is, using shared randomness. The local resources must 
be supplemented by non-local ones. A simple non-local resource is communication 
of information via classical bits (cbits), which we can quantify and thus provides a 
'measure of non-locality'. In this scenario, Alice and Bob try and output a and f3 
respectively, through a classical protocol, with the same probability distribution as 
if they shared the bipartite entangled quantum system and each measured his or 
her part of the system according to a given random Von Neumann measurement. 
^ ! As we have mentioned above, such a protocol must involve communication between 
Alice and Bob, who generally share finite or infinite number of random variables. 



The amount of communication is quantified [2f either as the average number of cbits 
C(P) over the directions along which the spin components are measured (average 
or expected communication) or the worst case communication, which is the max- 
imum amount of communication C W (P) exchanged between Alice and Bob in any 
particular execution of the protocol. The third method is asymptotic communica- 
tion i.e., the limit lirrin^ooC '(P n ) where P n is the probability distribution obtained 
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when n runs of the protocol carried out in parallel i.e., when the parties receive n 
inputs and produce n outputs in one go. A fundamental result for this scenario is 
that k2 n (k a constant) cbits of classical communication is required to simulate the 
correlations implied by a n qubit maximally entangled state [3|. This was followed 
by a remarkable result due to Toner and Bacon |4j who showed that a single cbit 
of communication is enough (apart from shared random variables) to simulate the 
correlations of a two qubit singlet state. Hence the nonlocalitv of two spin- 1/2 par- 
ticles in the singlet state is one cbit. We have later shown j^lo] that for simulating 
entanglement of arbitrary spin S singlet state, communication of n = log2{s + 1) 
cbits is enough, provided we confine only to spin measurements. 
Another resource which Alice and Bob can use to reproduce the quantum correla- 
tions is postselection. Here, Alice and Bob are allowed to produce a special outcome 
of their measurement. This corresponds to the physical situation where Alice and 
Bob's detectors are partially inefficient and sometimes do not click. Gisin and Gisin 
|7], inspired by Steiner's communication protocol, gave a protocol which simulates 
quantum correlations with shared randomness and with a probability 1/3 of abort- 
ing for either party. 

Another fruitful approach to this problem is to use PR-box jsj]. PR-box is a con- 
ceptual and mathematical tool developed to study non-locality, first proposed by 
Popescu and Rohrlich [9]. It was demonstrated in [10], that the correlations of the 
two qubit singlet can be simulated by supplementing hidden variables ( shared ran- 
domness) with a single use of the PR-box. Although mathematically based on the 
Toner and Bacon [4J result, this work is a major conceptual improvement, as the 
PR-box is a strictly weaker resource than a bit of communication, because it does 
not allow signaling. 

Finally, J. Degorre et.al Jl3|, [l4j have shown that the problem of reproducing the 
quantum correlation, in the worst case scenario for the two qubit singlet state and 
in the average scenario for the simulation of traceless binary observables on any 
bipartite state, with different resources, can be reduced to a distributed sampling 
problem. They have shown that the problem of reproducing the quantum corre- 
lations with different resources (classical communication, PR box, post selection) 
can be reduced to the problem of Alice and Bob agreeing on a sample from a dis- 
tribution depending on Alice's input. They introduced a method to carry out this 
distributed sampling in two steps. The first is the completely local problem of how 
Alice can sample a biased distribution depending on her input with the help of a 
shared uniform random source, and the second step is how Alice can share this 
biased sample with Bob by using communication, post-selection or non-local (PR) 
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box (for the two qubit state). 

It was soon discovered that for the general case of bipartite qubits in a partially en- 
tangled state two bits of communication is enough [4]. Brunner, Gisin and Scarani 



111 ] showed that a single use of PR-box is provably not sufficient to simulate some 
partially entangled two qubit states. One may think that the nonlocality of a par- 
tially entangled state shouldn't be larger than that of maximally entangled state. A 
reason for this apparently surprising result is that PR-boxes have random marginals, 
((a) = = (/?), a fact which is consistent with two qubit singlet state) while the cor- 
relations arising from partially entangled quantum states have nontrivial marginals. 
Thus it appears that it is especially difficult to simulate at the same time nonlocal 
correlations and nontrivial marginals, like these corresponding to partially entan- 



gled quantum states. Recently, Brunner, Gisin, Popescu and Scarani fl2| (hereafter 
referred to as BGPS) have given a procedure to simulate entanglement in non- 
maximally entangled states using four PR-boxes and one M-box. The M-box is 
defined by a b = [x < y] where [x < y] is the truth value of the predicate x < y for 
given values of x and y. In order to overcome the above difficulty due to nontrivial 
marginals, they introduce the concept of correlated local flips, which is independent 
of whether we use nonlocal boxes or classical communication to simulate entangle- 
ment. The idea is that some nonlocal box, or a classical communication protocol, 
first simulates non-local correlations with trivial marginals and then use local flips 
to bias the marginals. 

In this paper, we present a protocol to simulate the quantum correlations implied 
by a partially entangled two qubit state, in the worst case scenario, by using shared 
random variables as unit vectors uniformly and independently distributed over the 
unit sphere in K 4 and single use of a M-box and a PR-box. The paper is organized as 
follows. In Sec II, we recall the BGPS protocol. In Sec III, we present our protocol. 
Finally we conclude in section IV. 

II. BGPS PROTOCOL 

We briefly review the BGPS protocol and the idea of the correlated local flips. The 
problem is to simulate the quantum correlation implied by a general non-maximally 
entangled two-qubit state 

\^) = cos( 7 )|00) + sin( 7 )|ll) (0 < 7 < |). (1) 

Alice and Bob perform measurements along directions a and b on their qubits. Let 
a, (3 denote their outputs respectively. For binary outcomes (a, {3 G { — 1, +1}), the 



correlations are conveniently written as 

P(a, p\a, b) = i(l + aM A (a) + PM B (b) + a/3C(a, &)) (2) 

where 

M A (a) = ^aP{a,P\a,b) 

a,/3 

M B {b) = ^pP{a i P%b) (3) 

are the local marginals, and 

C(a, 6) = a(3P{a, P\a, b) (4) 

is the correlation term. Here we shall focus on pure non-maximally entangled states 
of two qubits. Thus the quantum correlation Pqm{&, /3\a,b) is given by 

M A (a) = ca z 
M B (b) = cb z 

C(a, b) = a z b z + s(a x b x - a y b y ), (5) 

where c = cos 27 and s = sin 27. 
This gives 

Pqm(®, P\a, b) = j(l + aca z + (5cb z + a/3C(a, b)) (6) 

where C(a, 6) is defined in Eq.([5J). In order to simulate this joint probability, BGPS 
12j proceed as follows. They first set up a procedure ( involving non-local boxes or 
communication) to simulate the joint probability Po(aj3\a,b) = |[1 + a/3Co\ where 
Co is the correlation (a/3)o = s ^2 l af3P{a^ f3) before flipping. Next, they invoked 
the local flip operation as follows. Alice (Bob) flips the output -1 with probability 
fa(fb) while the output +1 is left untouched. After the local flipping operation with 
probabilities f a and /& respectively by Alice and Bob, assuming /& > f a , the joint 
probability Po becomes 

P f (a, P\a, S) = |[1 + af a + 0f b + ap(f a + (1 - f b )C )]. (7) 

In order that Pf coincides with Pqm, they identify f a = ca z and ft, = cb z . Note 
that the condition f b > f a now becomes b z > a z . This gives, Co = a.B where 
B = (sb x , sb y , b z - c)/(l - c6j. If / > / 6j (a z > b z ) , then 

P f (a, p\a, S) = i[l + afa + Ph + afitfb + (1 - /o)C )]. (8) 
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where Co = A.b and A = (sa x , sa y , a z — c)/(l — ca z ). It is easy to see that A and 
5 are unit vectors. The joint probability P$ to be simulated before local flipping 
operation depends on whether (b z > a z ) or (a z > b z ) (via Co). But Alice (Bob) 
cannot have any information on b z (a z ). In order to pave way through this situation 
BGPS invoke a nonlocal box called M-box which has two real inputs x, y G [0, 1] 
and binary outputs m,n G {0, 1}. The M-box is defined by m n = [x < y] where 
[x < y] is the truth value of the predicate x < y for given values of x and y. 



III. THE PROTOCOL 



In this section, we present a protocol for simulating entanglement in an arbitrary 
non maximally entangled two qubit quantum state (Eq.flT])) in the worst case 
scenario, which makes a single use of M-box [12j and a single use of nonlocal box. 
The non-local box is used to obtain a random variable whose distribution depends 



on Alice's input |13|, \\A 



Let Alice and Bob share a PR-box and a M-box as well as the normalized vectors 
Ai, A2, A3, A<i, As, ^0 and Ai independently and uniformly distributed over the unit 
sphere in R 4 . The protocol runs as follows: 



(i) Alice and Bob input a z and b z values respectively in M-box, which outputs m 
and n (m, n G {0, 1}) as above. We replace m and n by p = 2m — 1 and q = 2n — 1 
so that j9, q G { — 1, 1}. Alice (Bob) gets p(q) without knowing q(p). 

(ii) Alice (Bob) prepares the unit vector u G M 4 ( v G M 4 ) by using shared 
random vectors jii, \x<i, /^3, ^4? As uniformly and independently distributed over the 
unit sphere in R 4 and p(q). These are given by 

u = -^~ u i H — u 2 (9) 

v = — - — vi H — v 2 (10) 



where 



ui,2 = (^1,2,^01,02) (11) 
vi,2 = (^1,2,^01,02), (12) 
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where i — \/—T and uok (^ofc) is the zeroth components of vector Uk 

Uip = sgn(c.jdi i 4 : )a + sgn{c.ji2$)A (13) 

uok = sgn{c.jl b ) \\\ u k \\ 2 - 1| 3 A; = 1,2 

#1,2 = sgn(c.fa i2 )b + sgn(c.p,i4)B (14) 



^ofc = |||^|| 2 -1|5 A: = 1,2. 



Here sgn(x) = 1 for x > and sgn(x) = — 1 for x < (x G M). 
i 6 I 3 and B E K 3 are defined in the previous section and c G 
vector shared between Alice and Bob. 



The unit vectors 
1R 4 is a fixed unit 



(iii) Now, Alice and Bob use two uniform shared random variables A G M 4 and 

A 

Ai G M. defined above and make a single use of the non-local (PR) box to output 

a = sgn(u.X s ) (15) 
ft = sgn(v.X s ). (16) 

A A 

where A s G R is a random unit vector distributed according to a biased distribution 
with probability density 



/ \ \ \u.X s \ 
Pu{X s ) = 



2tt 2 



13, 14 



(17) 
Note that 



and is not a shared random variable between Alice and Bob 
the distribution pa(X s ) depends on Alice's input via u. The corresponding protocol 
consists of Alice producing A s distributed according to Pu(X s ) (Eq.(0)) and then 
Alice and Bob make a single use of the non-local (PR) box so as to enable Bob 
to produce his output as in Eq.([T5j) without knowing A s . The protocol to do this, 
making a single use of non-local (PR) box, is given by a trivial modification in the 
'non-local box' protocol in Ref |13|, [l5 



We now find all the relevant averages at this stage of the protocol. Integrating 



over A s , we obtain: 

r 

< a(3 >\ = pu(X s )sgn(u.X s )sgn(v.X s )dX i 



1 



2tt 2 
1 

1 



u.Xslsgn^.X^sgn^.X^dX, 
u.X s )sgn(v.X s )dX s 
u. I X s sgn(v.X s )dX s 



2tt 2 

= u.v, (18) 



where the the first equality defines < a (3 >^ , §3 stands for the unit sphere in 



d4 



and the final integral is given by J g X s sgn(v.X s )dX s = (2tt 2 )v ( see lemma 4 in |l4j] ) 
Also it is easy to check that < a > x = =< (3 >~ x . 

The average before the local flipping operation (see below) (denoted (-)o), is ob- 
tained by averaging over {/}} ( denoted (-){^}). We get: 

< a > =< (3 >o= 



and 



< a/3 > =< u.v > {m} 

1+^1 + q „ A 1 -pi - q a a 
= — g ~ < >M +~ ~ < W 2-^2 >{ M } 

1 + pl — q _ 1 -pi + q _ 
= — g Y~ < Wl '^ 2 >M ~2 2~ < W2 ' Vl >{/i} 



(19) 



From Eq([TTj)-Eq([Tlj), we can calculate all Uip-Vip- For example: 

U1.V1 = U\.Vi — woi^oi = sgn{c.jli)sgn{c.il^)a.b 
+ sgn(c.jj2)sgn(c.jj,i)A.B 
+ sgn(c.fi2)sgn(c.£i3)A.b 
+ oJ 

- ^n(c./i5)|||wi|| 2 -l|^|||^i|| 2 -l|^ (20) 
Now by using < sgn(c.jli)sgn(c.fij) >{ M }= Sij and < sgn(c.fii) >{ fi }= 0, we obtain: 

< ui.vi > {fl} = < u 2 .v 2 >{ M }= a.B (21) 

< U1.V2 >{^} = < u 2 .vi >{ At }= A6. (22) 
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Substituting these in Eq(JT2D, we get: 

< a p > 0= l±Ela.B + l -=MA!b . (23) 

(iv) Alice and Bob perform the local flip operation with probabilities f a = ca z and 
fb = cb z respectively. We denote the averages after flipping described in section 
II by (•)/ and averages in the quantum state \ifj) (Eq.([T])) by (-)qm- Now, if fb > 
f a (b z > a z ) then p = q which implies from Eq.(EBj) that < a(3 >o= Co = d.5, 
so that < a/3 >/= ca z + (1 — cb z )a.B =< a/3 >qm- If f a > fb {a>z > b z ) then 
p q which implies from Eq.(|23j) that < a/3 >o= Cq = A. b, so that < a/3 > /= 
cb z + (1 — ca z )A.b =< a/3 >qm- The local flipping operation ensures that < a >/= 
f a = ca z =< a >qm and < j3 >/= fb — cb z =< j3 >qm- We see that the protocol 
simulates the joint probability Pqm{&, (3\a, b) as in Eq.(Ej). 

IV. CONCLUSION 

We have constructed a protocol to simulate quantum correlation implied by a 
partially entangled two qubit state, making single use of M-box and nonlocal (PR) 
box. The non-local box is used to achieve the required (local) outputs by the 
two parties which depend on a random variable not shared by the parties. This 
protocol simulates partially entangled two qubit state correlations in the worst case 
scenario. Our protocol is an improvement on previous protocols which used cbit 
communication (for example 2 cbit communication in Toner and Bacon model |4J) 



and non-local boxes (four PR-box and one M-box in BGPS model |12|), as it uses 
only one M-box and one PR-box which is a weaker resource. Obviously, this is a 
minimal resource in the worst case scenario. 
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